function [dat]=r_pp
% ---- Read PP parameter
% (mgC.m-2.d-1)

% ---- Create the Structure ...        
dat = struct('dummy',1);

dirPP='/calcul/pc-charria/OBSERVATIONS/SATELLITE/NPP/';

% In this case, EPPLEY method ... possibility to do with VGPM or CBPM
load([dirPP,'Satellite_NPP_1998_2000.mat'],'lon','lat','npp_eppley')
npp_eppley=permute(npp_eppley,[3 2 1]);

index_lon=find((lon <= 52) & (lon >= -104));
index_lat=find((lat <= 82) & (lat >= -21) );

dat.londata=lon(index_lon);
dat.latdata=lat(index_lat);
dat.data=npp_eppley(index_lon,index_lat,:);

n=1;
for y=1998:2000
    for m=1:12
        dat.timedata(n)=datenum(y,m,15);
        n=n+1;
    end
end

% ---- Clean the Structure ...
dat=rmfield(dat,'dummy');    
